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Abstract 

We study the nature of electronic states in one-dimensional con- 
tinuous models with weak correlated disorder. Using a perturbative 
approach, we compute the inverse localisation length (Lyapunov ex- 
ponent) up to terms proportional to the fourth power of the potential; 
this makes possible to analyse the delocalisation transition which takes 
place when the disorder exhibits specific long-range correlations. We 
find that the transition consists in a change of the Lyapunov expo- 
nent, which switches from a quadratic to a quartic dependence on the 
strength of the disorder. Within the framework of the fourth-order ap- 
proximation, we also discuss the different localisation properties which 
distinguish Gaussian from non-Gaussian random potentials. 

Pacs numbers: 73.20.Jc, 72.15.Rn, 72.10.Bg 

1 Introduction 

In recent years, the interest for one-dimensional (ID) models with corre- 
lated disorder has been steadily growing, as it has become progressively clear 
that correlations of the random potential can deeply affect the structure of 
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the electronic eigenstates and endow ID disordered models with far richer 
transport properties than it was previously thought. 

For a long time it was believed that ID systems could not display com- 
plex features like the metal-insulator transition which takes place in three- 
dimensional (3D) models, since it was known that all one-electron states 
are localised in ID systems with totally random potentials, regardless of the 
strength of disorder (see |J and references therein). Further research, how- 
ever, led to the study of some variants of the standard ID Anderson model 
which exhibited a discrete set of extended states (see, e.g., @). All these 
systems were characterised by random potentials with short-range spatial 
correlations, a feature absent from the original Anderson model where the 
site energies are totally uncorrelated. 

Eventually, the role of long-range correlations was also investigated and 
it was shown that this kind of correlations can produce a continuum of ex- 
tended states ||. This finding was corroborated by the work of Izrailev 
and Krokhin who, using second-order perturbation theory, managed to es- 
tablish an analytical relation between localisation length and potential pair 
correlators, and used this result to show how specific long-range disorder cor- 
relations lead to the appearance of mobility edges in ID discrete models M. 
An experimental confirmation of this result was obtained by studying the 
transmission of microwaves in a single-mode waveguide with a random array 
of correlated scatterers ||. The results established in || and |4j] for discrete 
lattices were subsequently extended to continuous models and related to par- 
allel phenomena occurring in different fields like the propagation of waves in 
random media |J and the dynamics of classical stochastic oscillators [0. 

The discovery that ID disordered systems can display a metal- insulator 
transition analogous to the one which takes place in 3D models constituted a 
crucial advancement for the understanding of anomalous transport properties 
of ID random models. This theoretical progress is also relevant from the 
technological point of view, because it paves the way for the construction of 
ID devices with pre-determined mobility edges which can be used as window 
filters in electronic, acoustic, and photonic structures The importance of 
the result makes therefore highly desirable to reach a complete comprehension 
of the link between long-range correlations of the disorder and the appearance 
of a continuum of extended electronic states in ID models. 

At the analytical level, our understanding of this problem rests on the 
results which were first derived in [|J. Unfortunately, both this work and the 
ones that have followed in its wake suffered from two main limitations: the 
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disorder was supposed to be weak, so that a perturbative approach could be 
applied, and all the analytical results were obtained in the second-order ap- 
proximation. (These constraints were absent from the work [[J of de Moura 
and Lyra, but they drew their conclusions from numerical calculations and 
did not provide any analytical insight on the delocalisation mechanism.) An- 
alytical calculations are greatly simplified by truncating the expansion of the 
inverse localisation length to the second-order term; the main drawback of 
this choice is that it leaves open the question of the true nature of the states 
of the "extended" phase. In fact, a vanishing second-order inverse localisa- 
tion length can be related to various physical phenomena: it can indicate 
that the electronic states are completely delocalised, but it can also be a sign 
of weak forms of localisation, characterised by power-law decay of the elec- 
tronic probability distribution; finally, it is the result to be expected when the 
spatial ranges of electronic states in the "extended" and "localised" phases 
differ by orders of magnitude. 

The present work constitutes an attempt to partially remove the limits of 
the previous analyses; in this study we still focus our attention on the case of 
weak disorder but, with the help of a systematic perturbative technique, we 
manage to go beyond the second order approximation and obtain analytical 
results correct to the fourth order of perturbation theory. By this way we can 
shed light on the true nature of the states that are classified as "extended" in 
the framework of the second-order theory and we are also able to discuss the 
differences which emerge at this refined level of description between Gaussian 
and non-Gaussian random potentials. 

In the case of Gaussian disorder with long-range correlations, it turns 
out that the electronic states are exponentially localised on both sides of 
the mobility edge identified in [[|; however, the dependence of their inverse 
localisation length on the disorder strength changes from a quadratic to a 
quartic form upon crossing the critical threshold. For weak disorder, this 
implies that electrons in the "extended" states can actually move over far 
longer distances than the electrons which find themselves in a "localised" 
state. The use of terms like "mobility edge" and "delocalisation transition" 
is therefore legitimate, provided that one keeps in mind that the qualifica- 
tion of "extended" must be understood as "localised on a large spatial scale" 
rather than "completely delocalised" . The analysis of the localisation length 
in the fourth-order approximation also reveals that, upon increasing the en- 
ergy of the electrons, a second threshold appears, which separates the states 
whose spatial extension scales with the inverse of the fourth power of the 
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potential from states that are characterised by an even weaker localisation 
(whose exact form cannot be ascertained within the limits of the fourth-order 
approximation) . 

The distinction between Gaussian and non-Gaussian disorder cannot be 
discussed within the framework of the second-order approximation, because 
the second-order inverse localisation length only depends on the second mo- 
ment of the random potential, and differences between Gaussian and non- 
Gaussian distributions only show up in the higher moments. A systematic 
computation of the higher-order terms of the inverse localisation length, how- 
ever, establishes a connection between the n-th term of the expansion and 
the corresponding n-th moment of the potential: therefore, the fourth-order 
results obtained in this work allow us a first analysis of the differences be- 
tween Gaussian and non-Gaussian disorder. Lifting the Gaussian require- 
ment seems to strengthen the localisation of electronic states, because the 
fourth-order inverse localisation length is increased by a new term propor- 
tional to the fourth cumulant of the random potential (cumulants are specific 
combinations of the moments which vanish in the Gaussian case but are oth- 
erwise non-zero). Specifically, we analyse in detail the difference between a 
Gaussian potential with long-range binary correlations and a non-Gaussian 
potential with the same pair correlators but with a fourth-order cumulant 
having exponential decreasing form: we find that the non-zero cumulant pro- 
duces quartic localisation of all electronic states that lie beyond the second- 
order mobility edge, thus wiping out the additional threshold that appears 
in the Gaussian case. 

This paper is organised as follows: in Sec. H we give an exact formulation 
of the problem and we define the model under study; in Sec. |] we expose 
the perturbative method used and we give the general results for the second- 
and fourth-order terms of the localisation length. We devote Sec. [| to the 
case of random potential with short-range correlations; the problem of long- 
range correlation is then discussed in Sec. [5|. Finally, we summarise the main 
results and express our conclusions in Sec. ||. 

2 Formulation of the problem 
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2.1 Definition of the model 

We consider the ID disordered model defined by the stationary Schrodinger 
equation 

h 2 

tf'lx) + eU(xU(x) = Ei/>(x) (1) 

2m 

where ip(x) represents the state of a quantum particle ("electron") of mass 
m and energy E moving in a continuous random potential U(x). The dimen- 
sionless parameter e is introduced to measure the strength of the disorder; 
for the sake of simplicity, in the rest of this paper we will adopt a system of 
units in which h 2 /2m = 1. We will set the zero of the energy scale at the 
average value of the potential and we will assume that E > on this scale. 

To complete the definition of the system under study, the statistical prop- 
erties of the random potential must be precised. In the first place, we will 
assume that model ([!]) is spatially homogeneous in the mean; from the math- 
ematical point of view, this requirement implies that the moments of U(x) 
satisfy the invariance condition 

(U{x 1 + 5)U{x 2 + 5)... U{x n + 5)) = (U{x 1 )U{x 2 ) . . . U{x n )) (2) 

for every spatial translation 5. In Eq. (|2|), and throughout this paper, the 
angular brackets denote the average over disorder realisations, or ensemble 
average. Property (Q) implies that the n-th moment of the random potential 
can depend only on the relative positions of the points (x\, x 2 , ■ ■ ■ , x n ); this 
allows one to represent the n-point correlator as a function of n — 1 relative 
coordinates according to the identity 

Xn(xi,x 2 , ■ ■ ■ ,x n _i) = (U(x)U(x + x 1 )U(x + x 2 ) ...Uix + Xn^)). (3) 

Our second basic assumption on model (|l|) is that its average features do not 
depend on the sign of the random potential, i.e., that the physical properties 
of the model are invariant under the transformation U(x) — > —U(x). This is 
ensured by assuming that all odd moments of the potential vanish 

(U(x 1 )U(x 2 )...U(x 2n+1 )) = 0. 

Finally, for reasons of mathematical convenience, we will suppose that a finite 
length scale l c exists for the model (|I|) such that the statistical correlations 
between the values U(xi) and U(x 2 ) of the potential become negligible when 
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the distance separating the points X\ and x 2 exceeds l c . For the particular 
case of the two-point correlation function this assumption translates in the 
condition 

(U{x 1 )U{x 2 )) ^0 for |xi - z 2 | » Z c I (4) 

more generally, the existence of a correlation length l c implies that the n- 
point correlation functions (U(xi)U(x 2 ) ■ ■ ■ U(x n ) must satisfy the so-called 
"product property", meaning that, if the points be divided 

in two groups with \xi — Xj\ ^> l c for Xi and Xj belonging to different groups, 
the average (U(xi)U(x 2 ) ■ ■ ■ U{x n ) takes the value obtained by averaging the 
two different groups separately. To deal with the case of long-range corre- 
lations, we will first make use of the previous hypothesis and derive results 
valid for any finite l c ; when possible, we will then take the limit l c — ► oo. 

Having thus defined the statistical properties of the random potential, 
we can proceed to enunciate the most specific assumption about the model 
under investigation. In the following, we will restrict our considerations to 
the case of weak disorder, as defined by the relation 

e 2 (U(x)U(x')) < E 2 . (5) 

and by the analogous conditions on the higher moments of the random po- 
tential. Notice that, to ensure that condition (|5|) is valid independently of the 
specific form of the function U (x), the parameter e must satisfy the condition 
e < 1. 

Eq. flU), together with the specified statistical properties of the random 
potential and the assumption of weak disorder, completely defines the model 
under study. Our goal consists in determining the spatial behaviour of the 
solutions of Eq. QTJ) ; to this end it is necessary to give a precise definition of 
the key parameter known as localisation length. 

2.2 Localisation length 

To shed light on the localisation properties of model ([[]), we will use as 
an indicator of the localised or extended nature of the electronic states the 
quantity 

A = lim — ln(^ 2 (x)k 2 + ^ ,2 (x)), (6) 

where the parameter k is related to the electronic energy by the relation 
E = k 2 . In definition (||) the wavefunction ip(x) can be supposed to be real, 
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so that no need arises to consider the complex extension of the logarithm 
function. 

In this paper we will always refer to the quantity as the inverse locali- 
sation length or Lyapunov exponent, although this parameter is more usually 
defined through the expression 

A = lim — (In U 2 (x)k 2 + ij' 2 (x))}. (7) 

x^oo 2,X \ J 

The parameters ([]) and (|7|) differ in principle, but both can be used effectively 
to ascertain the localised or extended character of the electronic states: in 
fact, they belong to the same family of generalised Lyapunov exponents 

which were introduced long ago for the study of localisation in ID disordered 
models (see and references therein). In technical terms, Eq. (|7]) defines the 
standard Lyapunov exponent L , whereas the parameter (|jD is equal to the 
"generalised Lyapunov exponent of order two" divided by a factor two, A = 
L2/2. Our preference for the definition (|6]) of the inverse localisation length 
over the possible alternative (|7|) stems mainly from the fact that the first 
choice makes the study of model ([I]) more amenable to analytical treatment. 
The exact relation between the Lyapunov exponents (§) and (0) is not a 
completely solved problem; however, it is clear that the two affine indicators 
can be used equivalently to determine whether an eigenfunction of Eq. ([[]) is 
localised or not. Indeed, for delocalised states the alternative definitions (||) 
and (0) give the same result, since both parameters vanish when the solutions 
of Eq. (P are extended. As for exponentially localised states, for the case of 
interest here, i.e., that of weak disorder, one can prove that the Lyapunov 
exponents (y) and (0) coincide at least to the second order of perturbation 
theory 0. In Sec. ^ we will come back to this point to show that the identity 
seems to hold also beyond the second-order approximation. 

From the physical point of view, the difference between the alternative 
parameters (||) and (0) is perhaps best understood if the problem of electronic 
localisation in model is re-defined in terms of the energetic instability of 
a stochastic oscillator. In fact, the physical comprehension of the model ([I]) 
can be enhanced by observing that solving the stationary Schrodinger equa- 
tion ([|) is completely equivalent to studying the dynamics of the classical 
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stochastic oscillator defined by the Hamiltonian 



The Hamiltonian (^j) represents an oscillator with frequency k perturbed by 
a noise the equivalence of this system with model (P can be seen by 
considering that the dynamical equation of the oscillator 

q(t)+k 2 (l + eZ(t))q(t) = (9) 

coincides with Eq. ([!]) if one identifies the electronic wavefunction ip(x) with 
the orbit q(t) of the oscillator and if the noise £(£) is related to the random 
potential U (x) by the identity 

Z(x) = ~U(x)/E. (10) 

Underlying the mathematical equivalence of Eqs. (JTJ) and (Bp is the phys- 
ical parallelism between the localisation of the quantum states of model ([I]) 
and the energetic instability of the stochastic oscillator (||) [0, [T(|. The con- 



nection between the two phenomena emerges clearly if one writes the inverse 
localisation length (§) in terms of the dynamical variables of the random 
oscillator. The Lyapunov exponent (01) can then be expressed as 



lhniln(fc 2 (p 2 (t) + g 2 (t))) 



A 

At 

and interpreted as the growth rate of the average energy of the random 
oscillator @. A positive value of the parameter (|6]) can therefore be read 
both as a sign of exponential localisation of the eigenstates of model (|I|) and 
as an indication that the stochastic oscillator (||) is energetically unstable. 
As for parameter (|7|) , it can be rewritten in terms of the oscillator variables 

as 

~ r T , x(t + 6) , 
A = hm hm / In — — at. 



r ^oo5-o7o x(t) 

In this form, the parameter (0) reveals itself as the exponential divergence 
rate of nearby orbits, i.e., as the Lyapunov exponent for the trajectories of 
the oscillator (H) . We can conclude that the parameters © and (0) are both 
growth rates, measuring respectively the increase of the oscillator energy and 
that of the distance between initially nearby orbits. 
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3 Perturbative expansion and general results 



To compute the localisation length @, one needs to determine the asymp- 
totic behaviour of (ip 2 (x)) and (ip' 2 (x)). A convenient way to achieve this 
goal consists in deriving a deterministic differential equation for the average 
squares of the wavefunction and of its derivative. To this end, we introduce 
the scaled derivative 

4>(x) = ij)\x)/k 

and the vector 

/ tp 2 {x) 
u(x) = (j) 2 (x) 

\ tp(x)(f)(x) 

Taking the Schrodinger equation ([!]) as starting point, it is easy to prove that 
u(x) obeys the stochastic equation 



du . . J _. 

— = (A + e£(x)B)u, 

ax 



(11) 



where £(x) is the scaled potential ( [TOP and the symbols A and B stand for 
the matrices 



A 




and B 



/ 

V -k 




We now need to replace the stochastic equation flTT|) with an ordinary differ- 
ential equation for (u(x)); this can be done in the following way. For a fixed 
initial condition u(0) = u , the average solution of Eq. flTTp can be formally 
expressed as 

(u(x)) = Z(x)u (12) 

with 



Z(x) = exp (Ax) 



1 + (Texp r eM{x')dx') 
Jo 



where Texp indicates the T-ordered exponential and M(x) is the random 
matrix 

M(x) = £{x) exp (-Ax) B exp (Ax) . (13) 
Differentiating both sides of Eq. (|1^) one obtains 

d(u) 



dx 



Z'(x)u Q ; 



(14) 
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on the other hand, by inverting the relation (|T2|), one arrives at 

uo = Z" 1 (x) (u(x)) . 

Substituting this result in the right-hand side (r.h.s.) of Eq. (0), one finally 
obtains the desired equation for the average of w(x) 

= Z , (x)Z- 1 (x)(u(x)) = K(x)(u(x)). (15) 

ax 

For weak disorder (e — > 0), the generator K(x) that appears in Eq. fllTf ) 
can be represented by an expansion in powers of the perturbative parameter e 

oo 

K(x) = J>"K n (x). (16) 

n=0 

Van Kampen |TT| and Terwiel |12| have shown that all partial generators 
K n (x) can be expressed in terms of specific combinations of the moments 
of the random matrix (|H|) known as "ordered cumulants" and they have 
provided a systematic set of rules to construct them. Applying Van Kampen's 
prescriptions to the present case, the zeroth-order generator turns out to have 
the simple form 

K = A 

which corresponds to neglecting altogether the effects of the weak random 
potential. As for the first-order term, it vanishes like all odd-order generators 
K 2ri+ i(x) because of the condition that the odd moments of U(x) be null. 
The second-order generator is less trivial and represents the first term where 
the effects of the random potential appear 

PX 

K 2 (x) = exp (Ax) / dx 1 (M(x)M(x 1 )} exp (-Ax) . (17) 
Jo 

Straightforward passages allow to derive from this compact expression the 
matrix elements of K^x), whose explicit form is 

(K 2 ) n (x) = (K 2 ) 12 (x) = (K 2 ) 13 (x) = (K 2 ) 23 (x) = (K 2 ) 32 (x) = 

(K 2 ) 21 (x) : 



1 r x 

le Jo X2 ^ ^ + cos ^ 2ky ^ dy 
i r x 

p J Q X2{y) [-1 + cos {2ky)\ dy 
K 2 ) 31 (x) = ^J\ 2 (y) sin (2ky)dy (18) 



K 2 ) 22 (x) = (K 2 ) 33 (x) 
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where the symbol \2 represents the two-point correlator of the potential U (x) 
defined by Eq. (|3|) . The fourth-order generator has a more complex form and 
depends both on the four- and the two-point correlators of the potential 

K 4 (x) = Jq dxi Jq 1 dx 2 Jq 2 dx 3 exp (—Ax) 
x [{M(x)M(x 1 )M(x 2 )M(x 3 }) 

-<M(a;)M(a; 1 ))(M(ar2)M(x3)) (19) 

-(M(x)M(x 2 ))(M(x 1 )M(x 3 )) 

- (M(x)M(x 3 ))(M(x 1 )M(x 2 )}] exp (Ax) 



The integrand in Eq. (19) is a non-trivial example of ordered cumulant; notice 
that it differs from an ordinary cumulant because the matrices M(xi) and 
M(x2) do not commute for x\ 7^ x 2 . The explicit form of the matrix elements 
of K 4 (x) can be worked out from Eq. (0); the expressions, however, are quite 
lengthy and do not add much to the physical understanding of the problem, 
so that we will omit them here. In theory, the Van Kampen scheme allows one 
to compute the terms of the expansion (|1^) to any desired order; in practice, 
however, the involved calculations become rather cumbersome beyond the 
second order. In this paper we will not consider terms of higher order than 
the fourth; within this approximation, we will be able to derive the standard 
second-order localisation length and the first correction to this result. 

In principle, the generators K. n (x) are functions of the spatial coordi- 
nate x (and of the initial condition u(0) chosen to solve Eq. (|TT|) ), but the 
dependence on these factors dies out as soon as the condition x 3> l c is ful- 
filled ||11|| . To determine the asymptotic behaviour of (u(x)), therefore, we 



need not bother with the space-dependent generators K n (x) but we have 
only to consider their constant asymptotic limits 

K n = lim K n {x). (20) 

Notice that the existence of the asymptotic generators (|20f) is mathematically 
ensured by the sufficient condition that the random potential U(x) possess 
a finite correlation length l c . This implies that, when we relax this condition 
in order to study the case of long-range correlations, we will have to pay 
attention on whether the asymptotic generators ( p0|) continue to be well- 
defined in the limit l c —>■ 00. 

To sum up, on the basis of the previous considerations we can conclude 
that, in order to analyse with fourth-order accuracy the asymptotic behaviour 
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of (u(x)), we can legitimately replace in Eq. ([D^) the exact generator with 
its asymptotic and truncated form 



K(x) ~ A + e 2 K 2 + e 4 K 4 . (21) 

The eigenvalue of the matrix ([H]) with the largest real part determines the 
exponential growth rate of the vector (u(x)) and, consequently, the inverse 
localisation length (|6|). The original problem is thus reduced to that of solving 
the secular equation of the matrix (pl|) ; determining the largest root of this 
equation, one obtains that the desired inverse localisation length has the form 

A = e 2 X 2 + e 4 X 4 + o(e 4 ) 

and that the general expressions for the terms A2 and A4 are 

1 r°° 

A 2 (A;) = - 1T / X2(x) cos(2kx)dx (22) 
4fc z jo 



and 

VO = J^dx^ixi) jo°° ^2X2(^2) sm(2kx 2 ) 

+ ^ lim^oo J^ x) dx 1 dx 2 dx 3 {[x4^i,^i + x 2 ,x x + x 2 + x 3 ) 

- X2(xi)x2(x 3 )} [cos(2fcxi) cos(2A;x3) — cos (2kxi + 2kx 2 + 2^X3)] 

- X2(%i + x 2 )x2(x 2 + x 3 ) [cos (2kxi + 2kx 2 ) cos (2kx 2 + 2kx 3 ) 

- cos (2/cxi + 2kx 2 + 2kx 3 )) - ^2(^1 + x 2 + x 3 )x2(x 2 ) 

x [cos (2kxi + 2kx 2 + 2kx 3 ) cos {2kx 2 ) — cos {2kx\ + 2kx 2 )]} 

(23) 

where u(x) is the integration domain 



u i x ) — x 2, ^3) : < xi, < x%, < x 3 , xi + x 2 + x 3 < x} . 

Expression (p2|) reproduces the standard second-order formula for the 
Lyapunov exponent of model (|l|) (see, e.g., |]); the result shows that the 
second-order inverse localisation length for an eigenstate with wavevector 
k is proportional to the cosine transform of the binary potential correlator 
taken at twice the value of the wave vector. Not surprisingly, the fourth-order 
Lyapunov exponent (|23|) has a more complicated form, involving both sec- 
ond and fourth moments of the potential. The increasing complexity of the 
mathematical formulae reflects the different degree of sophistication of the 
second- and fourth-order approximation schemes. From the physical point 
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of view, this difference is perhaps better understood if one visualises locali- 
sation as an interference effect (see, e.g., jnj). In the second-order (or Born) 
approximation, localisation is seen as an interference phenomenon generated 
by double scatterings of the electron at two points separated by a distance 
over which the random potential has non-negligible correlations: hence the 
connection between the second-order localisation length and the two-point 
correlator established by formula (|22|) . The picture corresponding to the 
second-order approximation, therefore, includes an electronic wavefunction 
(with wavevector k) which is scattered first backward, generating a wave of 
intensity 0(e) and then is rescattered in the forward direction, with inten- 
sity 0(e 2 ). To go beyond the Born approximation, one must include the 
interference effects produced by multiple scatterings in the analysis of locali- 
sation. The fourth-order approximation, for instance, takes into account both 
quadruple scatterings and couples of double scatterings. In the first case, an 
electronic wavefunction is scattered back and forth four times, generating 
a forward-scattered wave of intensity 0(e 4 ) which interferes constructively 
with the incoming wave; in the case of double scatterings, one considers the 
interference between two twice-scattered waves, each of intensity 0(e 2 ). This 
is the reason why the term (|23| ) involves the four-point potential correlator 
as well as various combinations of double products of two-points correlators. 

To make the result for the fourth-order term ( [23] ) mathematically simpler 
and physically more transparent, it is opportune to express the fourth mo- 
ment of the potential in terms of the second moments and of the fourth-order 
cumulant, writing 

Xifa, X 2 , X 3 ) = X2(X 1 )X2(X 3 - X 2 ) + X2(X 2 )X2(X 3 - X X ) ^ 

+ X2(X 3 )X2(X2 - %i) + A 4 (x 1 ,x 2 ,x 3 ) 



where we have used the symbol A4 to represent the fourth cumulant. Eq. 
represents the generalisation for non-Gaussian disorder of the familiar iden- 
tity which allows one to express the fourth moment in terms of those of order 
two in the Gaussian case. In fact, in expression (03) the fourth moment is 



given as the sum of two parts: the Gaussian combination of the second mo- 
ments and the cumulant, which represents the non-Gaussian contribution 
and vanishes in the Gaussian case. This interpretation of Eq. ( p4|) makes 
quite natural to distinguish two components in the fourth-order Lyapunov 
exponent and to write 

X 4 (k) = \[ G \k) + \{ NG \k) (25) 
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where the first term in the r.h.s. of Eq. fl2"5|) represents the Gaussian part 
which depends only on the second moments of the random potential, whereas 
the second term is proportional to the cumulant and constitutes the spe- 
cific manifestation of the non-Gaussian nature of disorder on the localisation 
properties of model dH). Inserting (|24l) in formula (|23| ) one obtains that the 
explicit forms of the Gaussian and non-Gaussian terms are 

Xf\k) = jo°° ^1X2(^1) jo°° ^2X2(^2) sin(2A;a;2) 

-f^p lim^oo dx 1 dx 2 dx 3 {x2(^i + x 2 )x2(x 2 + x 3 ) 

x [cos(2A;xi) cos(2kx 3 ) — cos (2kx± + 2kx 2 ) cos (2kx 2 + 2kx 3 )] (26) 

+X2(x 2 )x2(xi +x 2 + x 3 ) [cos (2kxi) cos (2kx 3 ) 

+ cos (2kxi + 2kx 2 ) — cos (2kx± + 2kx 2 + 2kx 3 ) (1 + cos (2kx 2 ))}} 

and 

1 ( 

\i NG \k) = — Urn dxidx-idx^^xx, x l + x 2 , x x + x 2 + x 3 ) 

4fc 4 Ju(x) (27) 

x [cos(2/cxi) cos(2A;x3) — cos(2/cxi + 2fcx 2 + 2kx 3 )} . 

The Gaussian term ( f2~6"D can be further simplified with some calculus work, 
provided that the binary correlation function decays quickly enough at in- 
finity, which is ensured by condition (|j). In this way one can lower the 
dimensionality of the integral (^), reducing it to a much more manageable 
one-dimensional integral. Introducing the functions 



(f c (k,x) = / X2(y) cos(2ky)dy 

V , (28) 

<Ps(k,x) = / X2W sm(2ky)dy 



the final expression for the Gaussian part of the fourth-order Lyapunov ex- 
ponent can be written in the form 



Xf\k) = — <p a {k,Q)[<p e {k,Q)+4tp c {0,0)] 

1 r 00 
- [<p c {k, 0) + 2y? c (0, 0)] ^ X 2(x)x cos(2kx)dx 



dx 



+ -^ip s (k,0) J X2(x)xsin(2kx)dx 



(29) 



1 

4~F 



ip c (0, x)(p c (k, x) — y?c(0, x) cos(2fcx) 



As for the non- Gaussian term (p7|), the dimensionality of the integral can- 
not be lowered without formulating some specific hypothesis on the form of 
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the cumulant A 4 . A change of variable, however, allows one to cast expres- 
sion (^7|) in the form 

(NG) 1 f°° f X [°° 

A4 \k) = — dx dy J dzA 4 (y,x,z) ^ 
x [cos(2%) cos(2kz — 2kx) — cos(2kz)} 

which is more convenient for our later use. 

Formulae (^) and fl3H| ) constitute the central result of this section. They 
represent the general form of the Gaussian and non-Gaussian parts of the 
fourth-order Lyapunov exponent (|25|); by applying them to specific cases, we 
will be able to perform an accurate analysis of the localisation properties of 
model d). 



4 Disorder with short range correlations 

In this section, we analyse the localisation properties of random potentials 
characterised by two different kinds of short-range correlations: a delta- 
correlated potential (white noise) and a potential with exponentially decaying 
correlations. 



4.1 White noise 

A white noise is defined as a random potential whose cumulants have the 
form 

A n (x 1 ,x 2 ,...,x n _ 1 ) = a n 8{x x )5{x 2 ) . . .<J(z n _i) (31) 

(see, e.g., H]). In Eq. ([T!]) the a n are constants which determine the strength 
of the cumulants, while (xi, x 2 , ■ ■ ■ , £ n -i) represent relative coordinates of the 
points (x, x + Xx, x + x 2 ■ ■ ■ , x + x n -\). 

Let us consider the Gaussian case first. In this case, all cumulants vanish 
except the second, which coincides with the correlation function 

A 2 (x) = X 2(x) = a 2 5(x). (32) 



Substituting the binary correlator (p2|) in Eq. ( p2|) , one recovers the standard 
result for the second-order inverse localisation length 

Hk) = ^- (33) 
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To determine the fourth-order correction to this result, the first step consists 
in observing that in the case of white noise the functions (|28|) take the form 



<f c (x) 
ip a (x) 






for 


X 


> 


<x 2 /2 


for 


X 


= 


(To 


for 


X 


< 



(34) 



Substituting ( p2[ ) and (|3~4]) in expression fl2"9|), one arrives at the conclusion 
that the fourth-order Gaussian correction is zero 

\[ G \k) = 0. 

This peculiar result is a consequence of the delta-like nature of the corre- 
lation function. To understand this point, one needs to go back to Eq. flHf ) 
which shows that the partial generator K 4 (x) is an integral of the fourth 
ordered cumulant of the random matrix ([y|). As pointed out in Sec. ||, an 
ordered cumulant differs from an ordinary one because the matrices M(x) do 
not commute for different values of x. In the special case when the correla- 
tion function has the form (|32"D, however, the integrand in Eq. ( |i9|) vanishes 
unless all integration variables are equal to x. Consequently, the matrices 
in flT9| ) commute and the ordered cumulant reduces to an ordinary one, which 
is zero because of the Gaussian assumption. Actually, this reasoning can be 
extended to all generators K n (i) with n > 2 and one is thus led to the con- 
clusion that, in the special case of Gaussian white noise, the asymptotic limit 
of the generator ( |IED is exactly equal to 

K = A + e 2 K 2 . (35) 

The fact that K 4 (x) vanishes does not automatically imply that the fourth- 
order Lyapunov exponent must be zero; however it is possible to prove that 
this is actually the case by solving the secular equation for the exact gener- 
ator (|33|). In this way one arrives at the result 

8k 2 128A; 6 V ; 

which shows that no term of order 0(e 4 ) arises from the second-order gener- 
ator K.2 in the case of Gaussian white noise. 

If we drop the Gaussian assumption, the fourth-order generator no longer 
vanishes; however, assuming that the fourth cumulant has the form (pi]), one 
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finds that the non-Gaussian term (|3H| ) is zero. We are thus led to the conclu- 
sion that, in the case of white noise, all eigenstates of Eq. (JJ) are localised, 
because the second-order Lyapunov exponent ( |33"D is strictly positive for all 
values of the wavevector k. However, the specific nature of the correlations 
prevents quadruple scatterings of electrons from producing any net destruc- 
tive or constructive interference effect on the wavefunction, so that neither 
an enhancement nor a reduction of localisation show up in the fourth-order 
approximation. 



4.2 Exponentially decaying correlations 

We now turn our attention to the case of disorder with exponentially decaying 
correlations. Specifically, we assume that the two-point correlation function 
and the fourth-order cumulant have the form 



X2(x) = cr 2 exp (-(3 2 \x\) (36) 

and 

A A (x,y,z) = 04 exp [-/3 4 (\x\ + \y\ + \z\)) . (37) 

In Eqs. fl36|) and (0), the constants (3 2 l and (3± l represent the range of the 
correlation function and of the fourth cumulant, respectively; note that the 
two parameters may be different. 



Inserting the binary correlator (36) in Eq. (22), one obtains that the 
second-order Lyapunov exponent is 



This result confirms the general rule that random potentials with short-range 
correlations produce localisation of all the electronic eigenstates (with the 
possible exception of a discrete set of extended states). Note that the Lya- 
punov exponent ( |3"%| ) tends to zero in the limit (3 2 — > 0: physically, this means 
that the localisation of the electronic states becomes weaker and weaker as 
the range of the disorder correlations stretches over increasingly larger dis- 
tances. In the opposite case, i.e. when the range of the correlation function 
tends to zero, one can recover from Eq. (|3q) the result (|33|) already derived 
for the white-noise case. More precisely, one obtains the white-noise expres- 
sion if the limit (3% — > oo is taken while keeping constant the integral of the 
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second moment: f*™ X2{ x )dx = 6 2 . This condition implies that the correla- 
tion strength a 2 must scale as a 2 = d 2 (3 2 /2; when this relation is satisfied, the 
second-order Lyapunov exponent fl38|) reduces to the white-noise form ( p3| ) 
for p 2 00. 

The binary correlator (|36D also determines the form of the Gaussian part 
of the fourth-order Lyapunov exponent. For x > the functions fl28|) take 
the form 

2 

(f c {k, x) = — -2 [(3 2 cos(2kx) - 2/c sin(2/cx)] exp (-(3 2 x) 

4k + p 2 

<j 2 

<£ s (fc, x) = 2 [2k cos(2fcc) + /3 2 sin(2A;x)] exp (— (3 2 x) ; 



correspondingly, Eq. ( ^9|) becomes 



V ; 4 fc 4 (/? 2 2 + k 2 ) {(31 + 4k 2 f V ' 



This result coincides with the expression obtained in Ref. [T_5[ for the fourth- 
order term of the inverse localisation length (|7|); an important consequence of 
Eq. (|39|) , therefore, is that the identity of the localisation lengths (f|) and (0) 
is not restricted to the second order, but holds at least up to the fourth order 
when disorder is weak and Gaussian and the correlation function has the ex- 
ponential form fl5Ep. This coincidence further corroborates the conclusions of 
subsection ^]2| about the substantial equivalence of the generalised Lyapunov 
exponents (H) and (|7|). 

To complete the discussion of the fourth-order approximation, we have 
to compute the non-Gaussian term (|30|) when the fourth cumulant has the 
form fl3?p. A straightforward calculation leads to 

\( NG Uk) = — ^ (40) 

4 {) 8 k 2 ((3 2 + k 2 )((3 2 + 4k 2 )(9(3 2 + 4k 2 ) K } 



An important conclusion which can be drawn from Eq. fl40|) is that the non- 
Gaussian character of the disorder produces an enhancement of localisation. 
This increase may be negligible in the present case, because the second-order 
Lyapunov exponent (^) never vanishes; that is not the case, however, when 
long-range correlations of the potential come into play, as we will discuss in 
the next section. 
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Considered together, Eqs. (|59"D and (JO) show that the Gaussian and non- 
Gaussian parts of the fourth-order Lyapunov exponent are both decreasing 
functions of the wavevector k, with the same asymptotic behaviour for large 
values of k. This implies that the relative importance of the two terms is 
determined only by the relative magnitude of the cumulant strengths a 2 and 
<7 4 and of the decay constants (3 2 and /3 4 ; when the corresponding parameters 
are of of the same order of magnitude, neither component of the fourth-order 
Lyapunov exponent dominates the other. Comparing Eq. (|38| ) with fl39"D 
and fl40|) , one can also observe that for large values of k the second-order 
Lyapunov exponent decays as A2(A:) oc whereas the fourth-order cor- 

rection vanishes as A4 oc 1/k 8 . This shows how the Born approximation 
becomes more and more accurate as the energy of the electrons is increased, 
in agreement with the general principle that electrons of higher energy are 
less sensitive to the details of the random potential. 

Eqs. (^) and fl40|) also show that the both parts of the fourth-order Lya- 
punov exponent tend to vanish when the correlation lengths f3 2 and (3^ 1 
tend simultaneously to infinity, as one can expect considering that the local- 
isation effects become less pronounced as the range of disorder correlations 
is increased. Note, however, that the two localisation lengths need not be 
equal; in particular, it is possible to consider the case of a strongly correlated 
potential for which the correlation function is constant in space (i.e., (3 2 — 0) 
and the spatial variation of disorder only shows up through the fourth cu- 
mulant (^) with /3 4 > 0. In this case the terms (|3Bp and vanish and the 
inverse localisation length coincides with the non-Gaussian component fl40|). 

Finally, we observe that the two parts (|39|) and (|^) of the fourth-order 
Lyapunov exponent tend to zero in the white- noise limit, in agreement with 
the results of the previous subsection, as can be seen by taking the limits (3 2 — > 
00 and 04 — > 00 (with the constraints J X2{x)dx = a 2 and J A 4 (x 1; x 2 , x 3 )d 3 x = 



5 Disorder with long-range correlations 

In this section, we focus our attention on the case of disorder with long-range 
correlations; more precisely, we consider a potential with correlation function 
of the form 

. sin(2fc c x) , n , 

X2(x) = a 2 y c ' . (41) 

x 
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The correlator ( flip is substantially different from those considered in the pre- 
vious section, because it lacks any finite length scale beyond which it becomes 
negligible; as we will see, this feature strongly alters the localisation proper- 
ties of the electronic wavef unctions. To enhance the physical understanding 
of the problem, we will analyse first the physical properties which descend 



from the specific form (51) of the two-point correlation function, postponing 
the discussion of non-Gaussian effects until the end of this section. 



As a consequence of the slow decay of the function (41) , great attention 
must be paid when applying the formalism developed in Sec. || to the present 
case. In fact, in the derivation of the general results of Sec. |3] we avoided any 
mathematical inconsistency by conveniently assuming the existence of a finite 
correlation length; one could therefore question the chevalier application of 
the general formulae (22) and (|29|) to a case in which no correlation length 
can be properly individuated and condition (|) is not satisfied. To avoid 
these difficulties, we will study long-range correlators of the form fl41|) as the 
limit case for (3 —>■ + of the more general correlation function 



X2{x) 



sin(2fc c x) 
o" 2 exp 



x 



-I3x) 



(42) 



Since the correlator (|42| ) decays exponentially for (3 > 0, the results of Sec. y 
can be safely used for this model; the localisation properties for the particular 
case (|4lD can then be deduced by discussing the limit form of the Lyapunov 
exponent for (3 — > + . 

Having defined our approach, we can proceed to derive the second-order 
Lyapunov exponent for the correlation functions (|42| ) and (|41"|). To simplify 
the mathematical form of the results and get rid of absolute value signs, in 
the rest of this section we will suppose that k > 0. This is not a restrictive 
hypothesis, because the wavevector k enters model ([!]) only via the energy 
E = k 2 which is an even function of k; consequently, the Lyapunov exponent 
is also even in k. 

Inserting fl42|) in Eq. (f22"|), one obtains 



\ 2 ({3,k) 



Q-2 
8fc 2 



< arctan 



4/% c 



(3 2 + 4(Jfe 2 - k 2 ) 



n6 (4k 2 c - Ak 2 - (3 2 ) 1 (43) 



where 6(x) is the step function defined as 



9{x) 



if 
if 



x > 
x < 
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In the limit j3 — > + expression (|4*3| ) becomes 

. / 7 \ ! — ttt for < /c < k c / a a\ 

\ 2 {k) = { 8k 2 (44) 

[ for k c < k 

The behaviour of the second-order Lyapunov exponent ([|3]) and fl44]) is rep- 
resented in Fig. [3]. 
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Figure 1: Second-order Lyapunov exponent 
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Note that the limit (3 — > + entails a sharp qualitative change in both 
the correlation function and the second-order Lyapunov exponent. For the 
correlation function, the limit (3 — > + brings about a transition from an 
exponential to a power-law decay; correspondingly, the strictly positive Lya- 
punov exponent (f43|) is converted in the inverse localisation length (fP|) which 
vanishes when the wavevector k exceeds the critical value k c . Comparing ex- 
pressions ( fi3|) and (fHj), therefore, one can conclude that the passage from 
a disorder with finite-range correlations to a disorder with long-range corre- 
lations causes the emergence of a continuum of delocalised states, separated 
by a mobility edge from the localised ones 0. 

This is obviously an important conclusion, because it shows that a kind 
of Anderson transition can occur in ID disordered models; however, it is 
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not completely satisfactory, because it is valid only within the limits of the 
second-order approximation. More precisely, one cannot label as "extended" 
the states for k > k c without a certain margin of ambiguity, because the fact 
that the second-order term (fHf ) vanishes does not guarantee that the same 
be true for all higher-order terms of the Lyapunov exponent. It is therefore 
impossible to ascertain on the basis of Eq. (fP]) whether the electronic states 
for k > k c are really extended or, rather, localised in a different form or on 
a much larger spatial scale than the states for k < k c . To shed light on this 
point, it is useful to determine the fourth-order Lyapunov exponent (|29|) in 
the case of long-range correlations of the form (f|l]). Once more, we will first 
compute the fourth-order term (|29| ) for the correlator ( f42|) and then we will 
use this result to study the limit (3 — > + . 

When the two-point correlation function has the form (|42|) , the func- 
tions (EBf) become 



(To 

<p c (k,x) = -jIm{E 1 [l3x-i2(k + k c )x}-E 1 [Px-i2(k-k c )x]} 
(p 8 (k,x) = ^-Re{-E l [f3x-i2(k + k c )x\+E 1 [f3x-i2(k-k c )x}}, 



where Ei(z) is the exponential integral defined in the complex plane by 



Ei{z) 



oo e ~t 



dt with I arg z\ < n 



see, e.g., |16|). Correspondingly, the Gaussian part of the fourth-order Lya- 
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punov exponent takes the form 



Ai G) (/3, k) = ^-A- 



a 2 



+ 



Ak 4 I 16 
1 



P 



P 



P 2 + A(k - k c ) 2 (3 2 + A(k + k c ) 2 



71 



In 



k 

8 I p 2 + Ak 2 ^ p 2 + A(k + k c ) 2 



p 2 + A{k + k c f 
p 2 + A(k - k c ) 2 



+ 



+ 



p 2 + Ak 2 p 2 + A(k - k c ) 2 
k + k c 



+ 
+ 



P 

8 
1 

2(3 2 + A{k + k c ) 2 
1 k k c 

2p 2 + A{k-k c ) 2 
k 



In 



In 



P 2 + (fc + 2k c f 

p 2 + k 2 
p 2 + (k- 2k c f 



p 2 + k 2 



71 



ax(k) + a 3 (k) + -s 3 (k) 



7t 



ai(k) + a 4 (k) + —s 3 (k) + 7is 2 {k) 



7T 

2a\{k) - a 2 (&) + ^s 3 {k) - —s\\k) 



P 2 + Ak 2 

where the functions a-iik) are defined as 
a±(k) 
a 2 (k) 

a 3 (k) = arctan 
a^(k) = arctan 



1 Af3k c 

2 arctan p + 4 (jfc2 _ kj) 

1 4/3fc c 

— ctrCtctll 77; TTTT 

2 /9 2 + /c 2 - 4fc 2 



/3 2 + k 2 + 2kk c 
2PK 



P 2 + k 2 - 2kk c 



and the symbols Si(k) represent the step-functions 
si(k) = 



82(h) 
ss(k) 



1 if [3 2 + k 2 - Aki < 
otherwise 



1 if P 2 + k 2 - 2kk c < 

otherwise 

1 if p 2 + Aik 2 - k 2 ) < 
otherwise 
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Finally, the function Fp(k) which appears in the r.h.s. of Eq. fl4"o] ) is defined 
by the integral representation 



Fp{k) 



o\k jo 



X2O) bc(0, 0) - y? c (0, x)] sin(2kx)dx. 



For k > 2k c the function (|^) can be expressed through the series 



Fp{k) 



1 00 

— y 

2k tx n 2 



(3 2 + Ak 



f3 2 + 4k 



(3 2 + 4(Jfe - k c ) 2 



(3 2 + 4(k + k c ) 2 

cos(nS2) > sin(ri(5o 



cos(n5i) 



where 



So 



arctan 
arctan 
arctan 



2k c 

T 

2{k + k c 



(46) 



(47) 



For k < \jk 2 + (3 2 — k c , the function F(k) can be represented in the alternative 
form 

i?fh\ 1 f ( 2k °\* P 2 + 4(k + k c ) 2 

F p (k) = - arctan j In - - ^ _ ^ + G/J (*) 

where Gp{k) is the series 



1 00 f_iV* 





"(fc + fcc) 2 " 


n 


\k-k c f 




{ 


_ (3 2 + k 2 _ 




_ {3 2 + k 2 _ 





with coefficients 



00 -y 



sin 



(K 

(2n + /) arctan f 



(48) 



Note that, when /3 > y/Sk c , the convergence regions of the series (fiTD and fl48|) 



overlap, so that there is no need to use the integral representation (46) to 
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compute the values of the function Fp(k); the integral must be evalu- 
ated only if j3 < VSk c and for those values of k comprised in the interval 
y/k* + (3 2 -k c <k< 2k c . 

The behaviour of the fourth-order Lyapunov exponent is represented 
for k < 2k c in Figs. || and [| and for k > 2k c in Fig. |[ As shown by the 



Figure 2: Fourth-order Lyapunov exponent for k < 2k c and (3 > k c 
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figures, for large values of (3 the fourth-order Lyapunov exponent is a smooth 
and decreasing function of k, assuming small but positive values everywhere 
(except, of course, than in the low-energy limit k — > 0, where the results are 
not valid because the weak disorder condition (|) is no longer satisfied). Phys- 
ically this means that, when the correlation range of the disordered potential 
is short, the fourth-order Lyapunov exponent fl45|) is a negligible correction to 
the second-order term (f43D . The behaviour of X^\f3, k) becomes more com- 
plex for small values of j3. First, the function looses its monotonic behaviour 
and even starts to assume negative values in an interval comprised within 
the range < k < k c . Note that the fourth-order term can become nega- 
tive as long as the second-order term is positive; physically, this means that 
the fourth-order term reduces the localisation effect. Finally, for (3 -C k c , 
the function (/3, k) develops a pronounced negative minimum in a left 
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Figure 3: Fourth-order Lyapunov exponent for k < 2k c and (3 < k c 
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Figure 4: Fourth-order Lyapunov exponent for k > 2k c 
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neighbourhood of k c and a sharp positive peak in a right neighbourhood k c . 
For values of k close to k c , therefore, the fourth-order term can represent a 
conspicuous correction to the leading term when (3 <C k c \ since in this limit 



the second-order term tends to assume the discontinuous form (53), one can 
conclude that the effect of the fourth-order correction is that of smoothing 
the Lyapunov exponent in the region k ~ k c . 

The expression of the fourth-order Lyapunov exponent (|45|) is quite com- 
plicated, but it simplifies significantly in the case of greatest interest, i.e., in 
the limit (3 — > + , when it takes the form 



\f\k) 



4k 4 



71 



32k 



In 



k ~\~ k c 
k k r 







7r k 2 + kk c + k 2 c 
k c {k 2 - W) 
71 2k r — k 



8 k c (k - k c ) 



+ F (k) 
+ F (k) 



for < k < k c 



for k c < k < 2k c 



for 2k c < k 



(49) 



The function ( J4U| ) is plotted in Figs. ^] and ^. In Fig. the graph of A4 (/?, k) 



Figure 5: Fourth-order Lyapunov exponent for [3 = 
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Figure 6: Fourth-order Lyapunov exponent for (3 = and k > k c 
10 2 




with (3 = 10~ 2 k c is also plotted for comparison. As can be seen from Fig. || 
the fourth-order Lyapunov exponent (pE9| ) displays a qualitative behaviour 
close to that of the corresponding term (|45| ) for small values of (3. There are 
two main differences, though: first, the fourth-order Laypunov exponent ( f4"5| ) 
is strictly zero for k > 2k c and, second, it diverges for k — > k c . 

This divergence must be disregarded, however, because the limit (3 — * + 
cannot be legitimately taken in a neighbourhood of k = k c . To understand 
this point we recall that, as pointed out in Sec. the existence of the asymp- 
totic generators (|20| ) cannot be taken for granted when one drops the assump- 
tion that a finite correlation length exists. In the present problem of 
this kind arises for the second-order asymptotic generator (pf) which can- 
not be defined for k = k c in the limit (3 — ■> + . In fact, as can be deduced 
from formula (|I8|), the matrix element (K 2 )3i of the second-order asymptotic 
generator is proportional to the sine transform of the correlation function. 
The sine transform of the correlator fl42"D is a well-behaved function as long 
as (3 > 0; in the limit (3 — *> + , however, it diverges for k = k c . As a con- 
sequence, the general formulae (|29|) and (|30D can be applied to the study of 
disorder with long-range correlations fl4"T|) only for wavevectors k which are 
not close to the critical value k c . We can thus conclude that formula (^9|) 
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describes the correct behaviour of the fourth-order Lyapunov exponent ev- 
erywhere, except that in a small neighbourhood of k c . (Strictly speaking, 
the previous analysis also invalidates the second-order result ( pHf ) for k ~ k c ; 
this expression, however, descends from formula (0), which has a broader 
range of validity since it can be derived with alternative methods that do 
not require the existence of the asymptotic generator K 2 .) 

Having established the range of validity of the expression (f49"|) , we can now 
comment the physical meaning of the result. The first remarkable aspect of 
the fourth-order Lyapunov exponent (|49|) is that it takes positive values in 
the interval k c < k < 2k c where the second-order term ( f44| ) vanishes. This 
means that the inverse localisation length, which is a quadratic function 
of the potential for k < k c , assumes a quartic form when the wavevector 
exceeds the critical value k c . This feature clarifies the nature of the transition 
which occurs at k — k c : when the threshold is crossed, the electronic states 
continue to be exponentially localised, but over spatial scales which are much 
larger than those typical of the states with k < k c . The transition at k c , 
therefore, does not bring a change from exponential localisation to complete 
delocalisation, but an increase of order 0(l/e 2 ) in the spatial extension of 
the electronic states. Since we are considering the case of weak disorder, 
£ <C 1, this increase can be huge: the bold use of terms like "mobility edge" 
and "delocalisation transition" adopted in Ref. || is thus fully justified. 

A second relevant feature of the fourth-order Lyapunov exponent (f49"|) 
is that it vanishes for k > 2k c . This implies that, when the potential dis- 
plays long-range correlations of the form a second transition takes place 
at k = 2k c , with the electrons becoming even more weakly localised. The 
spatial extension of the electronic states beyond this second threshold can- 
not be estimated within the fourth-order approximation, which allows the 
only conclusion that for k > 2k c the inverse localisation length must be an 
infinitesimal of order 0(e 6 ) or higher. 

The emergence of a second threshold at k = k c is an effect of the Gaus- 
sian nature of the potential that we have been considering up to now. The 
situation is different for non-Gaussian potentials, because in this case lo- 
calisation is enhanced by the cumulant-generated term ([3(]). The specific 
form of the non-Gaussian term (^) depends on that of the fourth cumu- 
lant of the disorder; to discuss a concrete example, we will consider the case 
of a random potential with correlation function and fourth cumulant of the 
forms ([11]) and (|3TD, respectively. The fact that the correlation function has 
the form (41) implies that all moments of the potential decay slowly, so that 
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it is correct to speak of disorder with long-range correlations at all orders. 
On the other hand, the non-zero fourth cumulant fl37|) makes the fourth mo- 
ment of the potential different from what it would be in the Gaussian case. 
Since the cumulant (|37|) decays exponentially, the fourth moment of the po- 
tential significantly differs from its Gaussian form only in a limited spatial 
range (which can actually be very small if (3± is large); yet this difference, no 
matter how small, is enough to induce the localisation of all electronic eigen- 
states with k > 2k c with inverse localisation length equal to fl40|). We can 
therefore conclude that, in the case of disorder with long-range correlations 
of the form (pEip, non-Gaussian potentials generate a stronger localisation of 
the electrons that their Gaussian counterparts. 

To conclude this section, we remark all the results obtained so far for 
the localisation of electrons in a disordered solid can be transposed to any 
other model described by Eq. ([]]): in particular, the results of this section can 
be applied to the dynamics of the stochastic oscillator (|8|). In the language 
appropriate to this model, the results of this section can be expressed by 
saying that, when the stochastic perturbation of the frequency displays long- 
range temporal correlations of the form (f|l|), the behaviour of the oscillator 
can be dramatically altered by letting the unperturbed frequency k cross 
some critical thresholds. Specifically, the oscillator is energetically unstable, 
with energy growing with the rate (fH|), for frequencies lower than k c . As 
soon as the unperturbed frequency exceeds the value k c , however, the energy 
growth rate drops by a factor 0(e 2 ) so that the energetic instability of the 
oscillator manifests itself only on long time scales. Finally, for k > 2k c the 
oscillator becomes even more stable, with a further reduction of the rate of 
energy growth which may be at least of order 0(e 2 ) or of a simple factor 
0(e°) depending on whether the noise is Gaussian or not. 

6 Conclusions 

In this paper we studied the localisation properties of ID models with weak 
disorder, focusing our attention on the role played by spatial correlations of 
the random potential in shaping the structure of the electronic states. Us- 
ing a perturbative approach based on a cumulant expansion technique, we 
were able to derive both the standard second-order expression for the inverse 
localisation length and a new result for the fourth-order correction. The 
knowledge of the fourth-order term of the Lyapunov exponent makes pos- 
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sible to investigate the delocalisation transition which takes place when the 
disorder exhibits specific long-range correlations. The analysis of Sec. [| re- 
veals that this transition consists in a sharp change of the scaling law for the 
inverse localisation length: when the electronic wavevector reaches a critical 
value, the Lyapunov exponent switches from a quadratic to a quartic depen- 
dence on the disorder strength e and, correspondingly, the spatial extension 
of the electronic states increases by a factor 0(l/e 2 ). In the case of Gaussian 
potentials, a second critical value of the wavevector exists, beyond which 
the electronic states become even more delocalised, with Lyapunov exponent 
A = o(e A ) for e -> 0. 

This additional transition is absent in the case of non-Gaussian poten- 
tials, which cannot be distinguished from their Gaussian counterparts in the 
second-order approximation but whose specific features emerge when the 
description is brought to the refined fourth-order level. Specifically, non- 
Gaussian potentials show their effect through the addition to the Gaussian 
Lyapunov exponent of an extra term proportional to the fourth cumulant of 
the random potential. This cumulant-generated term strengthens the local- 
isation effect and maintains the quartic scaling of the Lyapunov exponent. 
Loosely speaking, the requirement that the disorder be Gaussian seems to 
act as a constraint which reduces the degree of disorder and therefore favours 
the delocalisation of the electronic states, in agreement with the general prin- 
ciple which links the degree of randomness with that of localisation of the 
electronic states. 

From this point of view, this paper can be seen as a further step towards 
the construction of a hierarchy of disorders according to the strength of the 
localisation effects they produce in ID models. In totally disordered systems, 
for which the potential has no spatial correlation, all electronic states are ex- 
ponentially localised; when the degree, if not the strength, of the disorder is 
somewhat reduced by the presence of short-range correlations of the poten- 
tial, a discrete set of extended states may appear; finally, when the potential 
exhibits specific long-range spatial correlations and thus the randomness of 
the model is further diminished, a continuum of states with large spatial 
extension emerges. This paper clarifies the features of the extended states 
which appear in the last case and shows that, in the class of disorders with 
long-range correlations, Gaussian potentials seem to be "less random" than 
non-Gaussian ones. 
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